function eq = find_agg(eq, param, glob, options)


q = eq.yf.*glob.sf(:,1).^(-param.eta)./param.C;

D_out = (sum(eq.L.*upsilonp(q,param).*q))^(-1); 

eq.D_out = D_out;

end